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ABSTRACT 

Time-resolved spectra during the cooling phase of thermonuclear X-ray bursts in low- 
mass X-ray binaries (LMXBs) can be used to measure the radii and masses of neutron 
stars. We analyzed ~ 300 bursts of the LMXB 4U 1636-53 using data from the Rossi 
X-ray Timing Explorer. We divided the bursts in three groups, photospheric radius 
expansion (PRE), hard non-PRE and soft non-PRE bursts, based on the properties 
of the bursts and the state of the source at the time of the burst. For the three types 
of bursts, we found that the average relation between the bolometric flux and the 
temperature during the cooling phase of the bursts is significantly different from the 
canonical F oc T 4 relation that is expected if the apparent emitting area on the surface 
of the neutron star remains constant as the flux decreases during the decay of the 
bursts. We also found that a single power law cannot fit the average flux-temperature 
relation for any of the three types of bursts, and that the flux-temperature relation for 
the three types of bursts is significantly different. Finally, for the three types of bursts, 
the temperature distribution at different flux levels during the decay of the bursts is 
significantly different. From the above we conclude that hard non-PRE bursts ignite 
in a hydrogen-rich atmosphere, whereas for soft non-PRE and PRE bursts the fuel 
is helium-rich. We further conclude that the metal abundance in the neutron star 
atmosphere decreases as the bursts decay, probably because the heavy elements sink 
faster in the atmosphere than H and He. 
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1 INTRODUCTION 

Thermonuclear, type-I, X-ray bursts are due to unstable 
burning of H and He on the surface of accreting neutron stars 
in low-mass X-ray binaries (LMXBs). During these bursts, 
the observed X-ray intensity first sharply increases, typically 
by a factor of ~ 10 in about 0.5—5 seconds, and after that 
it decrease s more or less exponentially within 10—100 sec- 
onds fe.g..lLewin et al.lll993l : IStrohmaver fc Bildstenll2006l : 
iGallowav et all 120081 ). The total energy emitted by an X- 
ray burst is typically ~ 10 39 ergs. Some X-ray bursts are 
strong enough to lift up the outer layers of the star. During 
these so-called photo spheric radius expansion (PRE) bursts 
(Basi nska et ah 1984), the radiation flux that emerges from 
the stellar surface is limited by the Eddington flux. The 
neutron-star surface origin is supported by the fact that 
the inferred emission area during a burst matches the ex- 
pected surface area of a neutron star, assuming that the 
thermonuclear flash expands to cover the entire star during 
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the radius expansion and cooling phases of the burst. The 
emission area can be estimated from the fitting of the energy 
spectra during the cooling tails of bursts. There are a num- 
ber of theoretical and obs ervational arguments that support 
this assumption (see, e.g. iFrvxell fc Wooslevl[l982l: iBildstenl 
ll995l : ISpitkovskv et afll2002l ; IStrohmaver fe Bildstedl20oli r 

One of the best studied sources of X-ray bursts is the 
LMXB 4U 1636-53. Also kno wn as V801 Ara, 4 U 1636-53 
was discovered with OSO-8 (|Swank et al.| [l976). and was 
subsequently studied in great detail using obser vations with 
SAS- 3, Hakucho, Tenma, and EXOSAT (see iLewin et ail 
1 1981 for a r eview). The orbital perio d of this binary sys- 
tem is 3.8 hr (|van Paradiis et ailll990l), and the spin per iod 
of the neutron star is 581 Hz llstrohmaver et al.l Il998al fbh . 
Using EXOSAT, iDamen et all j 19891 ) detected 60 bursts 
from this source between 1983 and 1986; from observa- 
tions with the Rossi X-ray Timing Explor er (RXTE), 172 
burst s were detected up to 2007 June 3 (|Gallowav et al.l 
2008), and mor e than 250 bursts including data taken af- 
ter that date l|Zhang et all |2009| ). Most of these X-ray 
bursts have standard, single-peaked, fast rising and expo- 
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nentially decaying light curves. Some single-peaked bursts 
show photospheric radius expansion, while a few other 
bursts from 4U 1636-53 show multi- peaked light curves, al- 
though these are not PRE bursts jlMaurer fc Watts] |200S| ; 



iBhattacharvva fc Strohmaverll2006al ; IZhang et alj|2009h 

In this paper, we compare the temperature and 
emission-area distributions in PRE non-PRE and double- 
peaked bursts during their cooling phase. 



2 OBSERVATION AND DATA ANALYSIS 

We analyzed all data available from the RXTE Proportional 
Counter Array (PCA) of 4U 1636-53 as of May 2010. The 
PCA consists of an array of five collimated proportional 
counter units (PCUs) operating in the 2—60 keV range. We 
produced 1-s light curves from the Standardl data (0.125- 
s time resolution with no energy resolution) and searched 
for X-ray bursts in the light curves. We used the 16-s time- 
resolution Standard-2 d ata to calculate X-r ay colours of the 
source, as described in IZhang et a l. (2009). Hard and soft 
colours are defined as the 9.7 — 16.0/6.0 — 9.7 keV and 
3.5 — 6.0/2.0 — 3.5 keV count rate ratios, respectively, and 
intensity as the 2.0 — 16.0 keV count rate. The colour-colour 
diagram (CD) of all observations of 4U 1636-53 is shown 
in Figure [T] The position of the source on the diagram is 
parameterized by t he length of the solid curve 5 a (see, e.g. 
|MindezeEaL|[l999). The 5 a length is normalized to the dis- 
tance between S a — 1 at the top right vertex, and S a = 2 at 
the botto m left vertex of the CD. Similar to the Sz length in 
Z sources (|Vrtilek et al-lfiglxT h 5* a is considered to be a func- 
H|h 



tion of mass accretion rate l|Hasinger fc van der Klis| [i989). 

In order to study the bursts in detail, we analyzed the 
PCA Event data, E_125us_64M_0_ls, in which each individ- 
ual photon is time tagged at a ~ 122 jj,s time resolution in 64 
energy channels between 2 — 60 keV. We used all PCUs that 
were operating at the time of the X-ray burst to produce 0.5- 
s resolution light curves in the full PCA band. For the time- 
resolved spectral analysis of the bursts we extracted spectra 
in 64 channels every 0.5 s from the Event data of all avail- 
able PCUs. We corrected each spectrum for dead time using 
the methods supplied by the RXTE team. We generated the 
instrument response matrix files for each observation, and 
we fitted the spectra using XSPEC version 12.4.0. Because 
of the short exposure, in this case the statistical errors dom- 
inate, and therefore we did not add any systematic error 
to the spectra. We restricted the spectral fits to the energy 
range 3.0 — 20.0 keV. For each burst we extracted the spec- 
trum from the persistent data just before (or after) the burst 
to use as background in our fits; this approach, used to ob- 
tain the net emission of a burst, is well establ ished as a stan- 
dard procedure in X-ray burst analysis (e.g. iKuulkers et al.l 
2002). We fitted the time-resolved net burst spectra with a 
single-temperature blackbody model (bbodyrad in XSPEC), 
as generally burst spectra are well fitted by a blackbody. 
During the fitting, we kept t he hydrogen column density 
iV H fixed at 0.36 x 10 22 cm" 2 (|Pandel et al.ll200Sl ). and to 
calculate the radius of the emitting area in km, we assumed 
a distance of 5.95 kpc (|Fiocchi et al.ll2006l ). The model pro- 
vides the blackbody colour temperature (Tbt>) and a normal- 
ization proportional to the square of the blackbody radius 



(Rbb) of the burst emission surface, and allows us to esti- 
mate the bolometric flux as a function of time. 

We note that this procedure fails if the blackbody emis- 
sion during the burst comes from the same source that pro- 
duces the blackbody emission seen in the persistent emis- 
sion, since the diffe rence between two blackbod y spectra 
is not a blackbody (|van Paradiis fc Lewinl Il986l ) . This ef- 
fect is significant only when the net burst emission is small, 
and therefore problems may arise only at the start and tail 
of the burst, when the net burst emission is comparable 
to the underlying pe rsistent emission (see the discussion in 
IKuulkers etal]|2002h . 

To test whether fitting the net spectrum of the bursts 
affects our results, we proceeded as follows: For a subsample 
of the bursts we fitted the persistent spectrum just before 
each burst with a blackbody and a power law. During the 
burst, we fixed the power-law index and normalization to 
the value of the persistent emission, and let the blackbody 
parameters to vary. Using this procedure we found results 
that are consistent with the ones we found using the stan- 
dard method therefore, in the rest of the paper, we used the 
results we obtained using the standard method 

We defined the cooling phase for different types of 
bursts observed in 4U 1636-53 as follows: The cooling phase 
of non-PRE single-peaked and PRE bursts starts from the 
moment that the fitted blackbody temperature starts to de- 
crease until the end of the burst, defined as the time when 
the burst flux is within 3 a of the pre-burst persistent emis- 
sion. For double-peaked bursts we used two different defini- 
tions of the cooling phase, one that starts from the moment 
the blackbody temperature starts to decrease after the first 
peak of the burst until the end of the burst, and the other 
one from the moment the blackbody temperature starts to 
decrease after the second peak of the burst until the end 
of the burs t. We did not cons ider the triple-peaked burst 
reported in lZhang et~ai1 i|2009T l 



3 RESULT 

3.1 Temperature distribution during the cooling 
phase of X-ray bursts 

In order to understand the cooling tails of PRE, single- 
peaked no-PRE, and double-peaked bursts, we fitted the 
time-resolved spectra with an absorbed blackbody and com- 
pared the fitting parameters. From the CD in Figure [T] and 
hardness-intensity diagram (HID) in Figure [5] we found 
that all the PRE bursts and double-pea ked bursts have a 
hard colour small er than 0.85 (see also iMuno et al. 1 12004 
IZhang et al1l2009D . We therefore divided the bursts in four 
groups based on their persistent hard colour and their time- 
resolved spectrum (see Figures Q] and [3J : (i) Hard non-PRE 
bursts; these are single-peaked non-PRE bursts that took 
place when the persistent hard colour of the source was 
larger than 0.85. (ii) Soft non-PRE bursts; these are single- 
peaked non-PRE bursts that took place when the persistent 
hard colour of the source was smaller than 0.85. (iii), PRE 
bursts; (iv) Double-peaked bursts. 

In Figure [3] we investigate the relation between the 
bolometric flux and the color temperature during the cooling 
phase of the bursts using 0.5 s bins (see section 2). In this 
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Figure 1. Colour-colour diagram of all RXTE observations of 
4U 1636-53 up to May 2010. The grey points represent the colours 
of the source from all available RXTE observations. Each point 
in this diagram represents 256 s of data. The black crosses rep- 
resent the colours of the persistent emission of the source at the 
onset of a single-peaked non-PRE X-ray burst. The filled blue 
circles indicate the same for PRE burst, and the red squares for 
double-peaked bursts. The position of the source on the diagram 
is parameterized by the length of the black solid curve 5 a . The 
horizontal line divides hard and soft state of the source (see text 
for details). The diagonal solid bar divides the PRE bursts into 
two groups (S a > 2.1 and 5a < 2.1; see text for details.) 
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Figure 2. Hardness-intensity diagram of all RXTE observations 
of 4U 1636-53. The symbols are the same as in Figure [T] The 
PRE bursts are divided into two groups (intensity > 0.12 Crab 
and intensity < 0.12 Crab) by the vertical solid bar. 



so-called flux-temperature (FT) diagram we also show lines 
of constant inferred radius, calculated from F/(aT ). The 
top panels (a) and (b) of Figure [3] show the cooling phase of 
non-PRE bursts and PRE bursts, respectively; the bottom 
panels (c) and (d) show the cooling phase of double-peaked 
bursts after the first and second peak, respectively. In or- 
der to compare different types of bursts at the same source 
state, in panel (a) of Figure [3] we only plot the soft non- 
PRE bursts. At the beginning of the cooling phase (upper 
left part of the panels), both PRE and soft non-PRE bursts 
appear to have similar colour temperature. However, at the 
end of the cooling phase (lower-right part of the panels) the 



soft non-PRE bursts appear to have a larger spread of tem- 
peratures than the PRE bursts. Comparing the two bottom 
panels in Figure [3] we notice that there are some points at 
high flux level but below the 8-km radius red line in panel 
(c) that are not present in panel (d). These points are the 
ones between the two peaks in double-peaked bursts. After 
the second peak, in double-peaked bursts most of the points 
appear above the red line. 

In order to study the distribution of temperatures in 
the cooling phase of the bursts quantitatively, we divided 
the data of Figure [3] in four different groups based on flux 
levels. Panels (a) to (d) in Figure [4] show the distribution 
of the colour temperature during the cooling phase of hard 
non-PRE bursts (black), soft non-PRE bursts (red), PRE 
bursts (green) and double-peaked bursts after the second 
peak (blue), at different flux levels. In Table [T] we give the 
average and standard deviation of the distribution of tem- 
peratures for each flux level. We find that, as the bolometric 
flux decreases in the cooling tail, the average temperature 
of the hard non-PRE bursts becomes significantly higher 
than that of the soft non-PRE bursts. We also find that, as 
the flux decreases, the average temperature of soft no-PRE 
bursts becomes significantly higher than that of PRE bursts 
and double-peaked bursts. At higher flux levels, the PRE 
bursts dominate the sample. 

We carried out a Kolmogorov-Smirnov (KS) test to as- 
sess whether any two distributions are consistent with being 
the same. The results of the KS tests are shown in Table 
[5] We find that, except at the highest flux level, it is un- 
likely that any two sample distributions in Figure [4] come 
from the same parent population. We note that there are 
only eleven double-peaked bursts in our sample and there- 
fore the KS-test results are less reliable when we compare 
other samples to this type of bursts. Since it is apparent 
that the average temperature of different types of bursts is 
different (Table [1} , we repeated the KS test after we first 
aligned all distributions such that they all have the same 
average temperature. The new KS-test results are shown in 
Table [3] As expected, after aligning the distributions, the 
probabilities of the KS test increased. But at low flux level, 
the probability that two samples come from the same parent 
population remains small. This indicates that in these cases, 
not only the average, but also the shape of the distributions 
is different. 



3.2 Flux-Temperature relation during the cooling 
phase of X-ray bursts 

The evolution of the bolometric flux and the blackbody tem- 
perature during the cooling tails of all X-ray bursts in our 
sample appear to follow a power-law relation. We fitted a 
power law, Fb = ocT^, to the data for PRE, soft non-PRE 
and hard non-PRE bursts, separately, where Fb is the bolo- 
metric flux, Tbb is the blackbody colour temperature, a is 
the normalization and 7 is the power-law index. 

For all type of burst, we first selected data with jp- > 2, 
where <5F\, is the error of the bolometric flux, and we rebined 
the data by a factor 10 (we also tried j^- > 1 and j^- > 
3, but the power-law fitting results are consistent with the 
ones we found using j^h > 2 ). We then divided the data 
into three groups: data with Fb > 1 x 10 -8 erg cm -2 s , 
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Figure 3. Bolometric flux as a function of blackbody colour temperature for different types of bursts in 4U 1636-53 during the cooling 
phase of the burst. Each point on the plot represent 0.5 seconds. The diagonal lines represent lines of constant inferred radius. Typical 
90% confidence level error bars are indicated in the lower panels.. 



Table 1. The average and standard deviation of the temperature distribution for different types of bursts at different flux levels. 



Flux range (10 8 erg cm 2 s 1 ) 






0-1.0 


1.0-2.0 


2.0-4.0 


4.0-7.0 


hard non-PRE 


< 


kT bb > (keV) 


1.42 


1.79 


2.05 


2.37 






a (keV) 


0.27 


0.15 


0.18 


0.64 


soft non-PRE 


< 


kT bb > (keV) 


1.21 


1.61 


2.01 


2.42 






a (keV) 


0.23 


0.21 


0.20 


0.38 


PRE 


< 


kT bb > (keV) 


1.09 


1.58 


1.98 


2.36 






cr (keV) 


0.23 


0.20 


0.19 


0.19 


double-peaked 


< 


kT bb > (keV) 


0.98 


1.46 


1.83 


2.45 






cr (keV) 


0.23 


0.23 


0.38 


1.20 



F b > 2 x 10 -8 erg cm' 2 s" 1 , and all data, and fitted each 
group separately. The best fitting results are shown in Table 
[4] and the fits are shown in Figure [5] 

A power law fits well the data with flux larger than 
1 x 10 -8 erg cm" 2 s _1 for all type of bursts, but a power 
law does not fit all the data (without flux selection) for any 
of the three types of bursts (see Table [4]). We also found 
that: (i) During the cooling phase, the X-ray bursts of 4U 
1636-53 do not follow the expected F b oc T bb relation; (ii) 
within the same type of bursts, the power-law index changes 



significantly as the flux increases; (iii) within the same flux 
selection, the power-law index of hard non-PRE bursts is sig- 
nificantly higher than that of PRE and soft non-PRE bursts 
(except at the highest flux level where the fitting errors are 
large due to the limited data). Given that L b = aT^ h , this 
means that for a fixed flux level, hard non-PRE bursts have 
higher T bb than PRE and soft non-PRE bursts as the neu- 
tron star cools down or, for a fixed Tbb during the cooling 
phase, the apparent emitting area of PRE and soft non-PRE 
bursts is larger than that of hard non-PRE bursts. 
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Figure 4. The distributions of colour temperatures during the cooling phase for hard non-PRE bursts (blue), soft non-PRE bursts (red), 
PRE bursts (black) and double-peaked bursts after the second peak (green), respectively. Each panel is for a different flux range, as 
indicated. 



Table 2. KS test probabilities for the distributions in Figure HI 



Flux 


range (10 8 erg cm 2 s 1 ) 


0-1.0 


1.0-2.0 


2.0-4.0 


4.0-7.0 




soft 


non-PRE vs hard non-PRE 








3.0 x 10" 8 


7.3 x 10" 


l 




soft non-PRE vs PRE 





2.1 x 10" 6 


1.5 x 10" 1 


0.4 x 10" 


l 




hard non-PRE vs PRE 








6.2 x 10" 9 


9.8 x 10" 


l 


soft 


non-PRE vs double-peaked 





4.7 x 10" 12 


2.9 x 10" 2 


6.1 x 10" 


l 




PRE vs double-peaked 


1.1 x 10" 7 


3.2 x 10" 5 


3.9 x 10~ 2 


3.5 x 10" 


l 



3.3 Photospheric radius expansion bursts 

The PRE bursts concentrate on a narrow range of hard 
colours in the CD and HID, but they span a large range 
of soft colours in the CD, and a large range of intensity in 
the HID. We therefore divided the PRE bursts in two groups 
according to the position of the source in the CD and HID 
at the time the bursts started, one group with S a < 2.1 
(low-Sa PRE bursts) and the other with Sa > 2.1 (high-S^ 
PRE bursts; see the solid bar in Figure [I}. We also divided 
the PRE bursts in two groups according to the intensity of 
the persistent emission of the source at the time the burst 
started (see Figure [2]): the faint PRE burst, with an inten- 



sity lower than 0.12 Crabs, and the bright PRE burst, with 
an intensity higher than 0.12 Crabs. In panel (a) of Figure 
[6]we show the distribution of colour temperature during the 
cooling phase of bright PRE (black) and faint PRE bursts 
(red), and for low-S a PRE (red) and high-Sa PRE bursts 
(black) in panel (b) of Figure [6] in both panel we selected 
only data with flux less than 1 x 10 -8 erg cm -2 s . The KS 
test for these two groups for all flux levels, is shown in Table 
[5] The KS test still gives small probabilities at the lowest 
flux level; at other flux levels the samples are consistent with 
coming from the same parent populations. 
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Table 3. KS test probabilities after we have aligned all the distributions to the same average temperature. 



Flux range(10 8 erg cm 2 s x ) 


0-1.0 




1.0-2.0 




2.0-4.0 


4.0-7.0 




soft non-PRE vs hard non-PRE 


2.27 X 10 


-7 


3.0 x 10 _ 


5 


1.1 x lO- 2 


9.3 x 10- 


1 


soft non-PRE vs PRE 


4.6 X 10' 


-5 


1.1 x 10" 


2 


1.8 x 10- 1 


8.4 x 10- 


1 


hard non-PRE vs PRE 


5.8 x 10- 


11 


1.6 x 10- 


2 


1.9 x 10- 3 


9.5 x 10- 


1 


soft non-PRE vs double-peaked 


7.8 x 10- 


-1 


5.9 x 10- 


1 


1.3 x lO" 2 


7.9 x 10- 


1 


PRE vs double-peaked 


1.8 x 10" 


-2 


7.8 x 10- 


1 


4.5 x lO" 2 


8.3 x 10- 


1 



Table 4. Power-law fit for different types of bursts in the flux-temperature diagram. In the table, a and 7 are the normalization and 
the index of the best-fitting power law. 



Flux range 
10" S erg cm" 2 s" 1 


all 


> 1.0 


> 2.0 






7 


X Z /d.o.f 




7 


X*/d.o.f 




7 


X Z /d.o.f 


PRE 


0.28±n.01 


3.35±0.03 


420/182 


0.31±0.01 


3.22±0.06 


120/103 


0.30±0.03 


3.24±0.11 


42/55 


soft non-PRE 


0.28±0.01 


3.23±0.04 


627/256 


0.35±0.02 


2.91±0.06 


117/147 


0.34±0.03 


2.96±0.13 


26/44 


hard non-PRE 


0.097±0.004 


4.48±0.06 


719/321 


0.18±0.01 


3.58±0.12 


139/133 


0.26±0.06 


3.09±0.31 


49/44 




(a) all flux 




kT bb (keV) 

(b) flux > 1 









'e 








s 












X2 

3 














2.5 


2 

kT bb (keV) 

(c) flux > 2 





(g) all flux (h) flux > 1 (i) flux > 2 

Figure 5. The flux- temperature diagram of PRE (top panels), soft non-PRE (middle panels) and hard non-PRE bursts (bottom panels). 
From left to right the panels show, respectively, all data, data where the flux was larger than 1 X 10 -8 erg cm -2 s - , and data where 
the flux was larger than 2 X 1Q— 8 erg cm -2 s - 1 . The best-fitting power law for each selection is shown with the red lines. 
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Table 5. KS test probabilities for low 5 a (< 2.1) and high 5a (> 2.1) PRE bursts, and bright (I > 0.12 Crab) and faint (7 < 0.12 Crab) 
PRE bursts. 



Flux range (10 — 8 erg cm~ 2 s _1 ) 


0-1.0 


1.0-2.0 


2.0-4.0 


4.0-7.0 


low S a vs high 5 a 


7.6 x 10" 14 


5.0 x 10~ 2 


1.3 x 10~ 2 


3.6 x 10- 1 


bright vs faint 


9.1 x 10" 15 


3.6 x 10- 1 


4.8 x 10" 1 


1.3 x 10" 2 




kTJkeV) kTJkeV) 
(a) 0.0-1. OxlO -8 erg cm -2 s" 1 (b) 0.0-1. OxlO -8 erg cm" 2 s" 1 

Figure 6. The temperature distribution during the cooling phase of PRE bursts with Flux < 1 x 10 — 8 erg cm~ 2 s — 1 . Panel (a) shows 
PRE bursts with persistent intensity > 0.12 Crab (black) and < 0.12 Crab (red); panel (b) shows the PRE bursts with the S a value > 
2.1 (black) and > 2.1 (red). 



4 DISCUSSION 

During their cooling phase, the X-ray bursts in 4U 1636- 
53 do not follow the bolometric flux-temperature relation 
Fb oc T bb , which would be expected if the apparent emit- 
ting area on the neutron star remained constant during this 
phase. The relation between bolometric flux and tempera- 
ture is significantly different between photospheric (PRE) 
and non-photospheric (non-PRE) bursts, as well as between 
non-PRE bursts that happen when the source is in the hard 
state (hard non-PRE bursts) and the soft state (soft non- 
PRE bursts). We also found that the temperature distri- 
bution at different flux levels during the cooling phase is 
significantly different between PRE, hard non-PRE and soft 
non-PRE bursts. This is consistent with the fact that differ- 
ent types of bursts have different power-law indices in the 
flux-temperature diagram. 

also 



This result ( see 



iBhattacharvva et al.l l20ld ; ISuleimanov et al _ 
at variance with the findings of IGiiver et al.l ~ 201fj ). who 



Gottwald of all Il986l : 
~ 201Ch 



studied three PRE bursts in another source, 4U 1820-30, 
and found that th ese PRE bur s ts fol low the F b oc T bb 
relation quite well. iGiiver et all i|2010T I, however, did not 
consider two other PRE bursts in this source because 
the cooling phase of these two bursts showed a complex 
behaviour, in which the emitting neutron-star area ap- 
peared not to be constant. Here we studied 65 PRE bursts 
and 223 non-PRE bursts (out of which 94 and 129 were, 
respectively, hard non-PRE and soft non-PRE bursts), and 
we did not discard any burst in our sample. While some 
of the individual bursts in 4U 1636-53 could lay close to 



the F b oc T bb relation (see Fig. [3]), several others deviate 
significantly from that relation, and on average the relation 
is significantly different from Fb oc Tb b for any type of 
burst in our sample. The discr epancy between our results 
and those of IGiiver et all (|2010l ) could be due to differences 
between the two sources: 4U 1820-30 is an ultracompact 
binary, with an orbital period of 685s, in which the neutron 
star accretes H e from the companion (jStella et al. 1987; 
ICummindliool ). whereas 4U 1636-53 has a longer orbital 
period (3.8 hours; see §1), and the accreted material is 
mostly hydr o gen. It is also possible that the results of 
IGiiver et alj 1 2010l ) are biased by the small number of 
bursts they used, three PRE bursts in their case compared 
to 65 in our case, or by the fact that they ignored two of 
the PRE bursts in 4U 1820-30 for their analysis. 

Time-resolved spectra in the cooling phase of ther- 
monuclear X-ray bursts can be used to measure the 
radii and masses of neutron stars. The net spectra of 
the thermonuclear X-ray bu rsts are usually well described 
by a blackbody spe ctrum (|Strohmaver fc Bildstenl 120061 ; 
iGallowav et al.l 20081 ). which allows us to obtain R bb , T bb , 
the blackbody radius and colour temperature, respectively, 
from which we can calculate the bolometric flux of the neu- 
tron star. The departure from the F b oc Tb b relation in the 
cooling phase of all types of bursts in 4U 1636-53 means that 
the blackbody radius, R bb , is not constant during the cooling 
phase of the bursts. This could be due to either changes in 
the emitting area of the neutron star during this phase, or 
to changes in the colour-correction factor, which accounts 
for hardening of the spectrum arising from electron scat- 
tering in the neutron-star atmosphere (|London et al.|[l986l ; 



8 G. Zhang et al. 



iMadei et al.l l2003) : f c = T bb /T cff = y/R^/Ra,, where 
is the neutron-star radius observed at infinity. 

If f c remains constant throughout the cooling phase of 
all X-ray bursts, the differences of i?bb should due to the 
changes of the emission area on the neutron-star surface. 
From Figure [5] and Table [4] it is apparent that the emission 
area in hard non-PRE bursts is smaller than in PRE bursts 
and soft non-PRE bursts, and that at the end of the cooling 
phase, the emitting area decreases faster in hard non-PRE 
bursts than in the other types of bursts. iGottwald et al.l 
|l986T ) also found that in EXO 0748-676 the average R hh in 
the cooli ng phase incre ases as the persistent flux increases 
(see also [Bildstcn 1995]). If for hard non-PRE bursts only a 
fraction of the surface of the neutron star is emitting during 
the cooling phase, and the size of the emitting area decreases 
as the burst decays, one would expect to see burst oscilla- 
tions in the cooling tail of hard non-PRE bursts. However, 
most bursts with oscillatio ns in 4U 1636-53 occur when the 
source is in the soft state ijMuno et al.ll2004l ). 

Alternatively, if R^, should remain constant throughout 
the cooling phase of all the bursts, variations in i?bb can only 
be due to changes of f c . We can write the colour-correction 
factor as: 



fc 



Ro 



R(l + z) 



F 



(1) 



where d is the distance to the source, R is the neutron- 
star radius, z is the gravitational redshift, and a is the 
Stefan-Boltzmann constant. In FigureQwe plot f c as a func- 
tion of the neutron-star luminosity in Eddington units. For 
this plot we assumed R = 9 km, (1 + z) = 1.35, d — 6.0 Kpc, 
and we used equation |T]) and the same data (Tbb and Rhb) 
as in section [321 To calculate the fluxes in Eddington units, 
Fb/Fvdd = ib/iEdd, for PRE and soft non-PRE bursts we 
used the average peak flux of PRE bursts as the Eddington 
flux, -FEdd- Assuming that PRE and soft non-PRE bursts are 
He bursts, and hard non-PRE bursts are H bursts, for the 
hard non-PRE bursts we took 59% of the Eddington flux 
that we used in PRE bursts (assuming a hydrogen mass 
fraction, X = 0.7). 

Figure [7] shows that, as the luminosity increases, the 
colour-correction factor of PRE and soft non-PRE bursts 
first decreases, then increases, and finally remains more 
or less constant. For hard non-PRE bursts the colour- 
correction factor first decreases, and then remains more or 
less constant as the luminosity increases. At the same lumi- 
nosity level, f c in PRE and soft non-PRE bursts is smaller 
than in hard non-PRE bursts, while f c is consistent with 
being the same in PRE and soft non-PRE bursts. The data 
in Figure agr ee qualitatively wi th the predictions of the 
model of ISuleimanov et al.l f2010) (see their Figure 2), if 
we assume that during the cooling phase of hard non-PRE 
bursts the neutron-star atmosphere is H rich, whereas for 
PRE and soft non-PRE bursts the atmosphere is He rich. 
At high luminosity levels, the data show that f c stays con- 
stant, and does not increase as quickly as predicted by the 
model. We speculate that during the early phases of the 
cooling phase of all bursts, there is an enhanced abundance 
of heavy elements in the neutron-star atmosphere that were 
produced during the thermonuclear flash, and this effect is 



L t/ L Edd 

Figure 7. Colour-correction factor as a function of luminosity in 
Eddington units. The different colour points correspond to differ- 
ent types of bursts, as indicated in the plot. 



stronger in the PRE and soft non-PRE bursts wher e dredge 
up of ashes is more effective ((Weinberg et al.ll2005t ): as the 
neutron-star cools down, the heavy elements settle down 
faster than H and He, and therefore the heavy-element abun- 
dance in the neutron-star atmosphere dec reases as the burst 
decay s. This is reported by the models of ISuleimanov et al.l 
(2010), in which f c decreases as the high-element abundance 
in the neutron-st ar atmosphere inc reases. (The same is true 
for the model of iGiiver et al.l 120101 1 It is then possible that 
the differ ence at high luminosity b etween the data and the 
models of ISuleimanov et al.l l|2010l ) is due to changes of the 
chemical composition of the neutron-star atmosphere during 
the early cooling phase of the bursts. 

We also found a significant difference in the the dis- 
tribution of the blackbody temperature during the de- 
cay of X-ray bursts in 4U 1636-53 in different types of 
bursts at different flux leve ls (see also lin't Zand et al.ll2009l ; 
iBhattacharvva et alj |2010T ) . At low flux level, the average 
Tbb is higher in hard non-PRE bursts than in PRE and 
soft n on-PRE bursts (Figure [4] see also Ivan der Klis et all 
1990). In our sample we have ~ 130 soft non-PRE bursts 
and ~ 90 hard non-PRE bursts, but in panel (a) of Fig- 
ure [4] non-PRE bursts have more data points than soft 
non-PRE bursts. (Remember that each point in that Fig- 
ure represents a measurement every 0.5 s in the cooling 
phase of the bursts.) This means that hard non-PRE bursts 
have longer decay times than PRE and soft non-PRE bursts. 
The difference of Tbb and burst decay time can also be ex- 
plained by different c hemical composition of the fuel layer 
i|Gallowav et al.ll2008T ): If PRE and soft non-PRE bursts ig- 
nite in a He-rich environment, most of the fuel is burnt in 
a very short timescale, and therefore these bursts show fast 
cooling at the decay phase. The hard non-PRE bursts ignite 
in a H-rich environment. Hydrogen burning proceeds more 
slowly, because it is limited by the /3-decay moderated by 
the weak force. During the cooling of hard non-PRE bursts, 
the Hydrogen left in the atmosphere keeps burning. Also, 
since the hard non-PRE bursts ignite in a H-rich environ- 
ment, they have larger / c than the PRE and s oft non-PRE 
bursts that ignite in a He-rich environment (|Madei et al.l 
12004 ISuleimanov et al.ll201fj| ; IGiiver et al.ll20ia T" therefore 
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the hard non-PRE bursts show higher average Tbb than the 
PRE and soft non-PRE bursts. 

We note that the differences in the distribution of Tbb 
are more significant at low flux levels. These differences can 
also be due in part to differences in the underlying emis- 
sion of the neutron star due to continued accretion during 
the burst. Van Paradijs & Lewin (f986) argued that the 
net burst spectrum is not a blackbody, and at the end of 
the burst the Tbb of the net burst spectrum is directly re- 
lated to the neutron-star surface temperature just before the 
burst. Therefore, differences in the average Tbb may reflect 
differences in the neutron-star surface temperature, depend- 
ing on the state of the source at the time the burst starts. 
However, there are significant differences in the distribution 
of Tbb at low flux levels between PRE and soft non-PRE 
bursts, whereas both types of bursts occur when 4U f 636-53 
is in the soft state, and therefore the persistent spectrum be- 
fore the bursts should be similar. Also, both types of bursts 
should ignite in a layer with similar chemical composition, 
which makes it difficult to explain the differences by differ- 
ences in the pre-burst persistent emission, f c or emission 
areas. 



5 CONCLUSIONS 

We studied the cooling phase of type-I X-ray bursts in the 
LMXB 4U f636-53 using all available RXTE data. We di- 
vided the bursts in three groups according to their prop- 
erties and the spectral state of the source at the time the 
burst started: Photospheric radius expansion (PRE), hard 
and soft non-PRE bursts, respectively. (Soft non-PRE and 
PRE bursts occurred in the soft state of the source.) We 
found that, during the cooling phase of the bursts: 

• For all types of bursts, the average bolometric-flux vs. 
temperature relation, Tb — Tbb, is significantly different than 
the Tb oc T^b relation that would be expected if the apparent 
emitting area on the neutron star remained constant during 
the decay of the bursts. 

• For all types of bursts, the average Tb — Tbb relation 
cannot be described by single power law over the whole flux 
range. 

• The average Tb — Tbb relation is significantly different 
for PRE, hard non-PRE and soft non-PRE bursts. 

• In line with the previous conclusions, the temperature 
distribution at different flux levels is significantly different 
for PRE, hard non-PRE and soft non-PRE bursts. 

These results imply that either the emitting area on the 
neutron-star surface or, most likely, the colour-correction 
factor changes during the cooling phase of X-ray bursts. 

We calculated the colour-correction factor separately 
fo r the three types of b ursts. Compared to the models 
of ISuleimanov et al.l l|201(t ). the dependence of the colour- 
correction factor with luminosity (in Eddington units) is 
consistent with a scenario in which the main source of fuel in 
hard non-PRE bursts is hydrogen, whereas for soft non-PRE 
and PRE bursts the main source of fuel is helium. 

Based on the dependence of the colour-correction fac- 
tor with luminosity at high luminosity, we suggest that at 
the beginning of the cooling phase of the bursts there is an 
enhanced metal abundance in the neutron star atmosphere, 



and that the relative metal abundance decreases as the burst 
flux decreases. 
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